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We report on a study of charmless B meson decays to three-kaon final states. The results are 
obtained with a 78.7 fb -1 data sample collected on the T(4S) resonance by the Belle detector 
operating at the KEKB asymmetric energy e + e~ collider. The branching fractions for B decays to 
three-body K + K + K~, K°K + K~, K S K S K + , and K S K S K S final states are presented. We also 
make a first attempt to perform an isospin analysis of the three-kaon final states. 

PACS numbers: 13.25.Hw, 14.40.Nd 



INTRODUCTION 



Studies of three-body B decays can significantly broaden the understanding of B meson decay mechanisms and 
provide additional possibilities for CP violation searches. Analysis of the Kirir and KKtt final states is presented in 
Refs. Jj] ||]. In this paper we report results on the study of charmless B meson decays to three-kaon final states. The 
analysis is based on a 78.7fb~ 1 data sample, which contains 85.0 million BB pairs, collected with the Belle detector 
operating at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider || with a center-of-mass energy at the T(45) 
resonance. All results reported here are preliminary. 



THE BELLE DETECTOR 



The Belle detector Q is a large-solid-angle magnetic spectrometer that consists of a three-layer silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC) for charged particle tracking and specific ionization measurement 
(dE/dx), an array of aerogel threshold Cerenkov counters (ACC), time-of- flight scintillation counters (TOF), and an 
array of 8736 CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux return located outside the coil is instrumented to detect Kl mesons 
and to identify muons (KLM). Electron identification is based on a combination of CDC dE/dx measurements, the 
response of the ACC, and the position, shape and energy deposition of the associated ECL shower. 

Charged hadron identification is accomplished by combining the responses of the ACC and the TOF with dE/dx 
measurements in the CDC into a single value using the likelihood method: 

C(h) = C ACC (h) x C TOF (h) x C CDC (h), 

where h stands for the hadron type (p, K, ir). Charged tracks are identified as protons, pions or kaons by imposing 
requirements on the likelihood ratios (PID): 

PID < p) = mrw-v pm{K) = mrwf FID(,r) = W) " w - 1 - pwm 
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At large momenta (>2.5 GeV/c) only the ACC and dE/dx are used since the TOF provides no significant separation 
of kaons and pions. We use a GEANT based Monte Carlo (MC) simulation to model the response of the detector and 
determine acceptance [||. 

EVENT SELECTION 

The selection criteria are similar to those used in the analysis of B decays to the Kirn and KKir final states [p|. 
Charged tracks are selected with a set of track quality requirements based on the average hit residual and on the 
distances of closest approach to the interaction point in the plane perpendicular to the beam and the plane containing 
the beam and the track. We also require that the transverse track momenta be greater than 0.1 GeV/c to reduce 
the low momentum combinatorial background. For charged kaon identification we impose a requirement on PID(if), 
which has 86% efficiency and a 7% fake rate from misidentified pions. Charged tracks that are positively identified 
as electrons or protons are excluded. Since the muon identification efficiency and fake rate vary significantly with the 
track momentum, we do not veto muons to avoid additional systematic errors. 

Neutral kaons are reconstructed via the decay chain K°(K°) — > K$ — > ir + ir~ . The invariant mass of the two pions 
is required to be in the range \M(tt + tt') - M K o\ < 12 MeV/c 2 . The displacement of the 7r + 7r vertex from the 
interaction point (IP) in the transverse (r-<j>) plane is required to be greater than 0.1 cm and less than 20 cm. The 
direction of the combined pion pair momentum in the r-<j> plane is required to be within 0.2 rad of the direction from 
the IP to the displaced vertex. 

We reconstruct B mesons in the K + K + K~ , K$K + K~ , K$KsK + and KgKgKs three-body final states. The 
inclusion of the charge conjugate mode is implied throughout this report. The candidate events are identified by their 
center-of-mass (cm.) energy difference, AE = (JV E^ — E^, and the beam constrained mass, M\ x = tJ~E~F— (J2iPi) 2 , 
where Eh = \/s/2 is the beam energy in the cm. frame, and pi and Ei are the cm. three-momenta and energies of 
the candidate B meson decay products. In the first stage, we select events with Mbc > 5.20 GeV/c 2 and —0.30 < 
AE < 0.50 GeV, which is a larger range than that used in our previous report Q. This allows for more detailed 
studies of the background. For subsequent analysis, we also define a signal region of |Mb c — Mb| < 9 MeV/c 2 and 
\AE\ < 0.04 GeV and a AE sideband region defined as 0.05 GeV < \AE\ < 0.15 GeV. 

To determine the signal yield, we use events with Mb c in the signal region and fit the AE distribution to the sum of 
a signal distribution and an empirical background. The AE signal shape is parameterized by the sum of two Gaussian 
functions with the same mean. The widths and the relative fractions of the two Gaussians are determined from a MC 
simulation. We find that the signal MC simulation gives, in general, a narrower AE width than data. To correct for 
this, we introduce a scale factor that is determined from the comparison of the AE widths for B + — > D°ir + events in 
MC and experimental data; this correction is about 10%. The background from qq continuum events is represented 
by a linear function. The AE shape of the BB background is determined from MC simulation, as described below. 



BACKGROUND SUPPRESSION 

To suppress the combinatorial background from e + e~ — > qq continuum events, we use a set of variables that 
characterize the event topology. We require | cos# t hr| < 0.80, where 8 t hr is the angle between the thrust axis of the 
B candidate and that of the rest of the event; the distribution of | cos0 t hr| is peaked near 1.0 for qq and is nearly flat 
for BB events. We also use a Fisher discriminant ||, J 7 , formed from nine variables of a "Virtual Calorimeter" [Q, 
the angle between the candidate thrust axis and beam axis, and the angle between the B candidate direction and 
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TABLE I: Summary of results for B meson decays to three-body charmless hadronic final states. The branching fractions 
and 90% confedence level (CL) upper limits (UL) are quoted in units 10 -6 . For the modes with one neutral kaon the quoted 
reconstruction efficiency includes the K° — > Ks — » ir + it~ branching fraction. For modes with more than one neutral kaon only 
the Ks — * 7r + 7r~ branching fraction is taken into account. 



Three-body 




Efficiency 


Signal Yield 


B (90% CL UL) 


Reference 


Results from 


mode 




(%) 


(events) 






Ref. 


K+K+R- 




23.5 ± 0.50 


565 ± 30 


33.0 ± 1.8 ± 3.2 


This Work 


35.3 ±3.7 ±4.5 


K°K+K~ 




7.20±0.f7 


149 ± 15 


29.3 ±3.4 ±4.1 






K S K S K+ 




6.78±0.f9 


66.5 ±9.3 


13.4 ± 1.9 ± 1.5 






K s KsKs 




3.98±0.f7 


12.21" 


4.3l^ ±0.75 






K + K-ir + 




13.8 ±0.31 


93.7 ±23.2 


9.3 ±2.3 (< 13) 


i 


< 12 


K°K+n- 




4.53 ±0.16 


26.8 ± 16.6 


8.4 ±5.2 (< 15) 


i 




D°n+, D° - 


■ K+it- 


28.8 ±0.57 


4000 ± 66 








D~tv + , D~ - 




8.24 ±0.14 


451 ± 25 









beam axis. We make a requirement on T that rejects 53% of the remaining qq background while retaining 89% of the 
signal. 

We also consider backgrounds that come from other B decays. We conventionally subdivide this background into 
two types. One type is the so called "generic" BB background that originates from the dominant b — ► c tree transition. 
The description of these decays is taken from an updated version of the CLEO group event generator || . We find 
that the dominant background from this source is due to B — > Dh decays, where h stands for a charged pion or kaon. 
To suppress this background, we reject events where any two-particle invariant mass is consistent within 15 MeV 
(2.5a) with D° — > K + K~, D a — > K~tt + or D + — » K°tt + . We also reject events with a K + K~ invariant mass that 
is consistent with Xco — * K + K~ (\M(K + K~) — M Xc0 \ < 50 MeV). The other potential source of background is 
rare charmless B decays that proceed via b — > s(d) penguins or b — > u tree transitions. Since these final states are 
not included in the main generator decay table, they are generated separately. We studied a large set of potentially 
dangerous two-, three-, and four-body final states. We do not find any rare charmless B decay mode that produces a 
significant background to the three-kaon final states. 

RESULTS OF THE ANALYSIS 

The AE distributions for all three-kaon final states are shown in Fig. [j], where data (points with errors) are shown 
along with the background expectation (hatched histograms). To extract the three-body signal yields we fit the AE 
distributions. The results of the fits are summarized in Table [l| The statistical significance of the B° — > KsK$Ks 
signal, in terms of the number of standard deviations is 4.3a\ It is calculated as y— 2 ln(£o/£max), where £ max and 
Cq denote the maximum likelihood with the nominal signal yield and with the signal yield fixed at zero, respectively. 
The significance of the signal in all other three-kaon final states exceeds 10ct. For convenience, some of the results 
obtained in Ref. [H] are also given in Table § 

To determine branching fractions, we normalize our results to the observed B + — > D°n + , D° — > K + ir~ and 
B° — > D~tt + , D~ — * K°ir~ signals. This reduces the systematic errors associated with the particle identification 
efficiency, charged track reconstruction efficiency, and the event shape variables requirements. We calculate the 
branching fraction for B meson decay to a particular final state / via the relation 
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FIG. 1: AE distributions for B — > KKK three-body final states. Points with error bars are data; the open histogram is the fit 
result; the hatched histogram is the background. The straight line shows the qq continuum background contribution. 



B(B - /) = x B(B - Dn)B(D -» Kn), 

where Nf and are the numbers of reconstructed signal events for the final state / and that for the Dir reference 
process, and £/ and Ed-k are the corresponding reconstruction efficiencies determined from MC. We use the recently 
updated results on B + — > D°ir + and B° — > D~n + branching fractions from the CLEO Collaboration The AE 
distributions for the reference processes B + — * D°tt + , D° — * K + tt~ and B° — * D~ir + , D~ — * K°ir~ are shown in 
Fig. ||. The results of the fits are summarized in Table |[ 

The results of the three-kaon branching fraction measurements are presented in Table |. To determine the recon- 
struction efficiencies for K + K + K~ and K°K + K~ final states, we use a simple model jjj that takes into account the 
non-uniform distribution of signal events over the Dalitz plot. The three-body signal in this model is parameterized 
by a (f)K intermediate state and a fxK state, where fx is a hypothetical wide scalar state. For the KsK$K + and 
KsKsKs final states, the reconstruction efficiencies are determined from MC simulated events that are generated 
with a uniform (phase space) distribution over the Dalitz plot. The dominant sources of systematic error are listed 
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FIG. 2: AE distributions for B+ -> D°ty+, D° -» K + 7v~ (left) and B° -> D~ n + , D~ -> K°-k~ (right) events. The open 
histogram is the fit result and the hatched histogram is the background. The straight line shows the qq continuum background 
contribution. 



TABLE II: List of systematic errors (in percent) for the B — » KKK branching fractions. 



Source 


K+K+R- 


K S K+K+ 


KsK s K+ 


KsKsKs 


B — > Dtt and D —> Kn branching fractions 


7.7 


11.0 


7.7 


11.0 


Efficiency non-uniformity over the Dalitz plot 


2.2 


3.6 






Signal parameterization 


2.1 


5.6 


4.8 


7.6 


Particle identification 


4.0 


4.0 


2.0 


2.0 


Ks reconstruction 






5.0 


10.0 


MC statistics 


2.1 


2.4 


2.8 


4.3 


Total 


9.4 


13.7 


10.9 


17.4 



in Table ||. We estimate the systematic uncertainty due to variations of reconstruction efficiency over the Dalitz 
plot by varying the relative fractions of quasi-two-body states. The uncertainty in the reconstruction efficiencies due 
to the limited MC statistics is also included in systematic error. The uncertainty due to the particle identification 
is estimated using pure samples of kaons and pions from D° — > K~ir + decays, where the D° flavor is tagged using 
D* + — ► D°7r + decays. The systematic error due to uncertainty in the Ks reconstruction efficiency is estimated 
from the study of D* + — > D°ir + , D° — > Kstt + tt~ decays. We estimate the uncertainty due to the signal shape 
parameterization by varying the parameters of the fitting function within their errors. 

The BaBar Collaboration has recently presented results on charmless three-body B decays The reported value, 
B(B + — > K + K + K~) = (34.7 ± 2.0 ± 1.8) x 10 -6 , is in agreement with the result presented here and our previously 
published measurement Qj. 

To examine possible quasi-two-body intermediate states in the observed B — > KKK signals, we analyze the two- 
kaon invariant mass spectra. The Dalitz plot for B + — > K + K + K~ candidate events in the M^ C -AE signal region is 
shown in Fig. ||. Since there are two same-charge kaons in this case, we distinguish the K + K~ combinations with 
smaller, M(K + K~) m - m: and larger, M(K + K~) meiK , invariant masses. We avoid double entries by forming the Dalitz 
plot as M 2 (K + K-) ma ^ versus M 2 (K + K~) min . 

The K + K~ invariant mass spectra for events from the B signal region are shown as open histograms in Figs, ^(a)- 
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FIG. 3: Dalitz plot for B + — > K + K + K candidates in the B signal region. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



M ( K+K ")min (GeV/c 2 ) M ( K+K ") ma x (GeV/c 2 ) 

FIG. 4: Two-particle invariant mass spectra for B + — > K + K + K~ candidates from the B signal region (open histograms) and 
for background events in the AE sidebands (hatched histograms), (a) M(K + K~) m in invariant mass spectrum. The inset in 
(a) shows the 0(1020) mass region in 2 MeV/c 2 bins, (b) M(K + K~) ma ,^ spectrum with M(K + i\~) min < 1.1 GeV/c 2 and (c) 
M(if+A"-) max with 1.1 GeV/c 2 < M(K+K-) min < 2.0 GeV/c 2 . 

^(c). The hatched histograms show the corresponding spectra for background events in the AE sidebands, normalized 
to the estimated number of background events from the AE fit. The M(K + K~) m i n spectrum, shown in Fig. f|(a), is 
characterized by a narrow peak at 1.02 GeV/c 2 corresponding to the 0(1020) meson and a broad structure around 
1.5 GeV/c 2 that is consistent with a scalar state. In contrast to the B + — > K + tt + it~ three-body decay we also 
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FIG. 5: Dalitz plot for B — » KsK K candidates in the B signal region. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

M(K + K") (GeV/c 2 ) MfKgK*) (GeV/c 2 ) 



FIG. 6: Two-particle invariant mass spectra for B° — > KsK + K candidates from the B signal region (open histograms) and 
for background events in the AE sidebands (hatched histograms), (a) K + K~ invariant mass spectrum. The inset in (a) shows 
the 0(1020) mass region in 2 MeV/c 2 bins, (b) M(7y S ^' ± ) spectrum with M{K+K~) < 1.1 GeV/c 2 and (c) M(K.sK ± ) with 
1.1 GeV/c 2 < M(K + K~ ) min < 2.0 GeV/c 2 . There are two entries per B candidate in (b) and (c). 



observe a strong "non-resonant" enhancement in the K + K + K final state that extends over the full M(K + K ),„;„ 
mass range in Fig. ||(a). To plot the M(K + K~) m£ix mass spectrum we subdivide the M(K + A'~) m j n mass region into 
two ranges: M(K+K-) min < 1.1 GeV/c 2 and 1.1 GeV/c 2 < M(K+K-) min < 2.0 GeV/c 2 . The M(K + R-) mayi mass 
spectra for these two regions are shown separately in Fig. ||(b) and Fig. ^(c), respectively. The prominent two-peak 
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FIG. 7: Dalitz plot and two-particle invariant mass spectra for B — > KsKsK candidates in the B signal region. There are 
two entries per B candidate in (b). 



structure apparent in Fig. [|(b) is due to the 100% <p meson polarization in the B + — ► cf)K + decay, which results from 
angular momentum conservation. 

The Dalitz plot for B° — > KgK + K^ candidate events in the M^ C -AE signal region is shown in Fig. |[ The KK 
invariant mass spectra for these events are shown as open histograms in Fig. ^, and the hatched histograms show 
the corresponding spectra for background events in the AE sidebands. The structure observed in the K + K~ mass 
spectrum is very similar to that observed in M(K + K~) m - m spectrum for the K + K + K~ final state (see Fig. [|): a 
prominent peak that corresponds to the <f> meson and a broad enhancement in the higher K + K~ mass region. 

The numbers of reconstructed B + — > K$KsK + and B° — > KsK$Ks signal events are significantly smaller because 
of the additional suppression due to the K° — ► K$ — * 7r + 7r~ branching fraction. A Dalitz plot for the B + — > K$KsK + 
candidate events in the M^, C -AE signal region is shown in Fig. ^, along with two-kaon invariant mass spectra. It is 
apparent in Fig. [?] that the KgKg invariant mass spectrum is quite similar to the K + K~ mass spectrum observed 
in the K$K + K~ final state. Except for the absence of the 4> meson, which cannot decay to KgKs, we observe a 
similar broad structure in the higher KsK$ mass region. The low mass region of the K$K + mass spectrum shown 
in Fig. 0(b) agrees with background and exhibits no prominent structures. 



DISCUSSION & CONCLUSION 



Charmless B meson decays have attracted a considerable amount of attention in recent years, primarily because 
they provide a possible way to extract weak phases. An important check of the Standard Model would be provided 
by measurements of the same CP-violating parameter in different weak interaction processes. A good example is 
the comparison of the measurement of the coefficient of the CP violating sm(Amdt) term in the time dependent 
analysis of neutral B meson decays. In B° — > (cc)K° decays (where (cc) denotes a charmonium state) this coefficient 
is sin(20i). Precise measurements of sin(20i) (also known as sin(2/3)) have recently been reported by the Belle and 
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BaBar experiments |1(| . The best known candidates for b — > s penguin dominated processes where this quantity can be 
measured independently are B° — > <fiK° and B° — > rj'K decays. However, these modes have small branching fractions 
of order 10~ 6 — 1CP 5 (including secondary branching fractions). Thus, very large numbers of B mesons are required 
to perform these measurements. This is especially true for the <f>K a final state. The large signal observed in the three- 
body B° — > KsK + K~ decay mode, where the <j>Ks two-body intermediate state gives a relatively small contribution, 
would significantly increase statistics if these events could be used. There are two possible complications: (1) Whilst 
the 4>Ks state has fixed CP, the CP-parity of the three-body K$K + K~ final state is not fixed. If the fractions of CP- 
even and CP-odd components are comparable, the K$K + K~ state will not be useful for a CP violation measurement; 
(2) Possible b — > u tree contributions may introduce an additional weak phase in the B — > K$K + K~ amplitude 
and cloud the interpretation of any observed CP violation. The b — ► u contribution in B° — > 4>Ks is expected to be 
negligible (since (j> is a pure ss state), but this is not necessarily the case for the three-body KsK + K~ final state. 
Here we discuss the possibility of the use of the three-body B° — > K$K + K~ decay mode for CP measurements. 

The decays of B mesons to three-body Khh final states can be described by b — > u tree- level spectator and b — > s(d)g 
one-loop penguin diagrams. Although b — » u VF-exchange and annihilation diagrams can also contribute to these final 
states, they are expected to be much smaller and we neglect them in the following discussion. 

B meson decays to final states with odd numbers of kaons (s-quarks) are expected to proceed dominantly via the 
b — > sg penguin transition since, for these cases, the b — > u tree contribution has an additional CKM suppression. In 
contrast, decays with two final state kaons proceed via the b —> u tree transition with no b — > sg penguin contribution. 
This allows us to estimate the b — > u tree contribution to final states with three kaons via the analysis of KKir final 
states. It is also important to note that the b — > sg penguin transition is an isospin conserving process, while the 
b — > u tree and b —> dg penguin transitions are isospin violating. 

This is illustrated for B + — » K + K + K~ decay in Fig. |§|. The diagram that corresponds to the main b — > s penguin 
contribution is shown in Fig. ||(a). The b — ► u tree contribution, Fig. ||(b), has an additional Cabibbo suppression due 
to the W + — > su vertex. The corresponding diagram without Cabibbo suppression (W + — > du) is shown in Fig. ||(d) 
and expected to be the dominant contributor to the K + K~ir + final state. A quantitative estimate of the size of the 
b — > u tree amplitude is provided by the ratio 

B{B+^K+K-.+ ) w ff K \ 2 2 , 



F =\aSF?~b(b+^k+k + k-) x {t) xtanPC ' 

where is the total amplitude for the B + — * K + K + K~ decay and A^Jf^ ^ s its 6 — *■ u tree contribution. The 

{f k/ fir) 2 factor, where f„ = 131 MeV and fx = 160 MeV are pion and kaon decay constants, respectively, takes into 
account the corrections for SU(3) breaking effects in the factorization approximation, and 0c is the Cabibbo angle 
(sin B c = 0.2205±0.0018) 0. Using the results for B+ -> K + K~Tr + and B+ -> K+K+R- branching fractions from 
Table |, we obtain F — 0.022±0.005. Similarly, for B° decays to KgK + K~ and KsK + n~ final states, respectively, we 
find F = 0.023 ±0.013 (< 0.037), where the second number is obtained using the upper limit for the B° — > KsK + n~ 
branching fraction. The small value of F indicates that we can neglect the b — > u tree contribution to the B — > KKK 
rates and perform an isospin analysis of the three-kaon final states. From an isospin decomposition of B mesons wave 
functions we obtain the following relations between three-kaon branching fractions, 

B{B° -> K°K+K-) =B(B+ -> K + K + K~) x — ; (1) 

T~B + 

B(B° -> K°K+K-) = B(B+ -> K+K°K°) x ISl, (2) 

T~B+ 

where the factor Tgo /tb+ takes into account the difference in total widths of charged and neutral B mesons. We use 
the first relation, Eq. [I], as a check of our assumption that the isospin violating contribution is small, and calculate 
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FIG. 8: Diagrams for B+ — * K + K + K decay: (a) b — > s penguin; (b) and (c) b — * u trees, and for B + — > K+K -r+ decay: 
(d) and (e) b — > u trees; (f) 6 — > d penguin. 



the ratio 

i? , y x Z*± = »K*K + K- x £™_ x ^ = Q 95 ± Q n ± 

o(-B+ — > K+K+K ) r B o N K+K+K - e K o K+K - t b o 

where we use signal yields (iV) and reconstruction efficiencies (e) from Table [| instead of branching fractions to reduce 
the systematic error, and t b +/t b o = 1.091 ± 0.023 ± 0.014 |l2j. The calculated value agrees with unity within its 
statistical error. 

Let us now consider the second isospin relation, Eq. |2|, in more detail. The B+ -» K+K°K° decay results in three 
different observable states: K+ K$Ks, K+K^Kl and K+ KgK B . The relative fractions of these states depends on 
the relative fractions of states with even and odd orbital momenta in the K°K° system. Bose statistics requires that 
the K°K° wave function be symmetric (and, therefore, CP even), independently of the relative orbital momentum, 
I, of the neutral kaons. As a result, a K°K° system with even orbital momenta can only decay to K$Ks or K^Kl 
final states (with equal fractions), while a K°K° system with odd orbital momenta can only decay to the K$Kl final 
state. Thus, the K+K°K° wave function can be written in the following form 



|JTW >= a \K + KsK s > + \K+K L K L > + >> 

V2 

where a and /3 are unknown coefficients constrained by a 2 + (3 2 = 1. 

In this experiment, we observe only the K + K$Ks component of the K + K°K° final state. Measuring the B + — > 
K+K+K~ and B+ — ► K+KgKg branching fractions and using the isospin relation, Eq. ||, with the wave function 
decomposition, Eq. [|, we can determine the parameter a 2 , 

2 = n B( B+ -> K+K S K S ) tbo_ = Nk+KsKs £k°k+k- ^> U] 

B(B° -> K°K+K~) r B+ N KsK+K - e K+KsKs t b+ ' {) 
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Here, the parameter a 2 characterizes the fraction of states with even orbital momenta in the K°K° system in the 
three-body K + K°K° final state. At the same time, due to isospin symmetry, a 2 also gives the fraction of states 
with even orbital momenta in the K + K~ system of the three-body K°K + K~ final state. Since the total angular 
momentum of the K°K + K~ system is zero, the orbital momentum of the K + K~ pair relative to the remaining 
neutral kaon, I', is equal to I. Thus, the CP-parity of the K$K + K~ three-body system is (—1)', and a 2 also gives 
the fraction of CP-even component of the three-body K$K + K~ final state. Using the information from Table |, 
we obtain: a 2 = 0.86 ± 0.15 ± 0.05. Note that the K$K + K~ three-body final state includes the (f>Ks state which 
is CP-odd. We remove B° — > <j)Kg events by requiring \M(K + K~) — M<p\ > 15 MeV; the number of remaining 
K$K + K~ events is 123 ± 14. The value of the a 2 for remaining events is: ct 2 on ^ = 1.04 ± 0.19 ± 0.06. Assuming 
isospin symmetry, we can also use the B + — » K + K + K~ final state instead of B° — ► K°K + K~ to determine a 2 . In 
this case Eq. [| becomes 

a2 = B(B+ -> K+K S K S ) = N K+KsKs e K+K+K - 

B(B+ -> K+K+K-) N K+K+K - e K+KsKs ' 1 1 

which gives a 2 = 0.82 ± 0.12 ± 0.06 and a 2 on . = 0.97 ± 0.15 ± 0.07. The two ways of computing a 2 and a 2 on , are in 
good agreement. This is evidence for the dominance of the CP-even component in the three-body non-</> KsK + K~ 
final state. 

In conclusion, we have measured branching fractions for charmless B mesons decays to the three-kaon K + K + K~ , 
K$K + K~ , K$KsK + , and KsK$Ks final states. The isospin analysis of the three-kaon final states presented here 
implies several conclusions that should be carefully examined further. The three-body KsK + K~ final state may be 
a good candidate for a CP violation measurement in b — + s penguin transitions. In this case a significant increase (by 
a factor of about four) in the statistics as compared to the B° — > <j)(K + K~)K$ final state, which is a well known 
mode for the measurement of CP violation in b — + s penguin dominated decays, is possible. 
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